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Abstract. The NUGA project is a high-resolution (0.5"-l") CO survey of low luminosity AGN 
including the full sequence of activity types (Seyferts, LINERs and transition objects). NUGA 
aims to systematically study the different mechanisms for gas fueling of AGNs in the Local 
Universe. In this paper we discuss the latest results of this recently completed survey, which 
now includes newly acquired subarcsec resolution observations for all targets of the sample. The 
large variety of circumnuclear disk morphologies found in NUGA galaxies (m = 1, m = 2 and 
stochastic instabilities) is a challenging result that urges the refinement of current dynamical 
models. In this paper we report on new results obtained in 4 study cases for NUGA: NGC 4826, 
NGC7217, NGC 4579 and NGC 6951. 



1. The NUGA project: searching for AGN feeding mechanisms 

While it is widely accepted that the onset of activity results from the feeding of a 
supermassive black hole by the gas reservoir from its host galaxy, there is no consen- 
sus on which mechanisms can remove virtually all of the gas angular momentum and 
drive the infall down to scales of tens of pc. Furthermore, it is unknown whether these 
mechanisms are at work only in active galaxies, or, alternatively, if the key difference 
between AGN and pure starburst or quiescent galaxies resides in the availability of gas to 
supply the central engine. The NUclei of GAlaxies-NUGA-project (Garcia-Burillo et al. 
2003a; Garcfa-Burillo et al. 2003b) is the first high-resolution (~0.5"-l") CO survey of 
12 low luminosity AGN (LLAGN) including the full sequence of activity types (Seyferts, 
LINERs and transition objects). Observations were carried out with the IRAM Plateau 
de Bure Interferometer (PdBI). Thus far, no interferometric survey had been focused 
on the study of AGN in the Local Universe (e.g., Heifer et al. 2003). High-resolution 
observations are paramount to achieve a sharp view of the distribution and kinematics 
of molecular gas down to the critical scales for AGN feeding (<100 pc). On these scales, 
secondary modes, embedded in kpc-scale perturbations (e.g., bars, spirals and tidal in- 
stabilities), are expected to take over. Furthermore, to get a complete picture of the AGN 
feeding issue NUGA relies on a multi-wavelength approach. Information on the stellar 
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Figure 1. Gravitational stellar torques on molecular gas evaluated in the 600 pc disk of the 
LINER NGC4826 (adapted from Garcfa-Burillo et al. 2003b). The strength of the i- Fourier 
components of the stellar potential, Qj, derived from HST-NICMOS data is shown in the upper 
left panel. The lower left panel shows the radial variation of the average torque per unit mass 
(t(r)). The radial extent of Arm I and Arm II, both indicated in the right panel (which overlays 
isovelocity contours onto the 12 CO(1-0) intensity map of NGC 4826), is also given. AGN feeding 
seems to be inefficient at present (see text). 

potentials and the star formation processes in NUGA targets is obtained from available 
HST and ground-based optical/NIR images. This multi-wavelength strategy allows us to 
test the efficiency of AGN feeding mechanisms by a direct determination of the stellar 
gravitational torques exerted on the gas and, also, by numerical simulations of the gas 
dynamics. 

The first images issued from NUGA reveal a wide variety of morphologies in the cir- 
cumnuclear molecular disks of AGN hosts (Garcfa-Burillo et al. 2003a; Garcfa-Burillo 
et al. 2003c). Various gravitational instabilities develop at different spatial scales, in- 
cluding 77i=2, 777=1 and stochastic perturbations. Some galaxies host several coexisting 
instabilities, while others host mainly one type of instability. The weak correlation be- 
tween activity type and nuclear morphology of the AGN host suggests that the AGN 
duty cycle in LLAGNs is shorter than the relevant time scale of the instabilities present 
in the maps. The relevant question to be addressed is to what extent the identified per- 
turbations favour AGN feeding or, eventually, if some of them can inhibit the process. As 
an illustration of the non-obvious quest for evidence of ongoing feeding in AGNs, we re- 
port below on recent results obtained in 4 study cases for the NUGA project: NGC 4826, 
NGC 7217, NGC 4579 and NGC 6951. 

2. AGN feeding in LINERs: NGC 4826 and NGC 7217 

The circumnuclear molecular gas disk of the LINER/transition object NGC 4826 shows 
the prevalence of m=l instabilities (see figure ^ adapted from Garcfa-Burillo et al. 
2003b). Gas kinematics are perturbed by streaming motions related to m=l modes. The 
non-circular motions associated with the inner 777=1 perturbations (lopsided instability 
and inner m = 1 spiral-Arm-I) fit the pattern expected for a trailing wave developed 
outside corotation ('fast' wave). A paradoxical consequence is that the inner 777=1 per- 
turbations would not favour AGN feeding. An independent confirmation that the AGN 
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Figure 2. We overlay 12 CO(2-l) contours on HST NIR color images (darker shades for redder 
colors) for two Seyferts in the NUGA sample observed at subarcsec resolution at the PdBI: 
NGC4579 (L1.9/Syl.9) (Garcfa-Burillo et a/.2004, in prep) and NGC6951 (Sy2) (Schinnerer 
et a£.2004, in prep), appearing in panels A and C, respectively. The locations of the AGN 
are given by star markers; the detections of ~1O 6 M0 molecular gas components near the AGN 
(denoted as G) are indicated by the arrows. The kinematics of the gas components very near the 
AGN are dissimilar in these two Seyferts: highly non-circular in NGC4579, and mostly circular 
in NGC 6951. This is illustrated in panels B and D showing the major axis p-v plots (oriented 
as shown in left panels) for NGC 4579 and NGC 6951, respectively. In both cases, NUGA maps 
show tantalizing evidence of what could be ongoing AGN feeding episodes. 



is not being generously fueled at present is found in the low values of the gravitational 
torques exerted by the stellar potential for R<530 pc (figure^. While the radial variation 
of the stellar torques seem to account qualitatively for the changing signature of stream- 
ing motions, the maximum value of the mean torque is exceedingly small: ^50(fcm/s) 2 . 
The combination of the stellar perturbations seems to make the gas lose its angular 
momentum rather inefficiently (see Garcia-BuriUo et al. 2003b for details). 

NUGA observations of the LINER NGC 7217 have been compared with numerical 
simulations in Combes et al. (2004) The emerging scenario in this LINER is apparently 
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in direct contrast with what is found for NGC 4826: in NGC 7217, molecular gas is piled 
up in an axisymmetric ring and kinematics is seen to be dominated by regular rotation. 
N-body simulations developed in Combes et al. (2004) served to monitor the formation 
and evolution of the gas ring. Results from these simulations show that gas inside the 
ring experiences positive gravity torques from the stellar oval perturbation which builds 
up the resonant ring near the ILR; hence, gas in the inner 700 pc is experiencing an 
outward flow. As for NGC 4826 (but for totally different reasons!) we find no evidence of 
significant ongoing fueling of the NGC 7217 nucleus. Only a GMC-like non self-gravitating 
unit coincident with the AGN seems to betray a past accretion episode. 



3. AGN feeding in Seyferts: NGC 4579 and NGC 6951 

The S1.9/L1.9 barred galaxy NGC 4579 has been mapped at ^0.5" resolution (see 
figure b adapted from Garcfa-Burillo et al. 2004, in prep). The molecular gas distri- 
bution in the inner 500pc mimics the gas response to a bar potential. The 12 CO(2-l) 
map shows 1.6xl0 8 M Q of molecular gas piled up in 2 leading lanes. However, the m — 2 
point-symmetry in the gas distribution breaks up in the inner 100 pc disk where lopsid- 
edness seems to take over. As illustrated by figure [2K the distribution of neutral gas 
close to the AGN consists of a 150 pc off-centered ringed disk which is also identified in 
the V-I color HST image of the galaxy. There is little molecular gas near the AGN: the 
closest molecular complex (at a radius r^80pc) has a mass of ^a few 10 5 M©. Modelling 
work (in progress) is key to interpreting the complex kinematics of molecular gas near the 
AGN, characterized by the presence of highly non-circular motions (figure Et). Results 
from numerical simulated orbits will shed light on the still controversial feeding issue in 
this Seyfert. 

NGC 6951 is a prototypical Seyfert 2 galaxy for which subarcsec resolution CO maps 
have been recently completed within the NUGA project (figure|2t,d adapted from Schin- 
nerer et al. 2004, in prep). Molecular gas distribution in the inner 700 pc shows two highly 
contrasted nuclear spiral arms containing a significant gas reservoir of 4x 10 8 M Q which is 
presently feeding an ongoing star forming episode, also identified by its prominent radio- 
continuum emission (Ho & Ulvestad 2001). The m = 2 gas instability is likely reflecting 
gas crowding along the X2 orbits of the stellar bar. As can be guessed from figure [2t, 
only a small amount of molecular gas seems to have succeeded in making its way down 
to the AGN. A compact molecular complex of ~a few 10 6 Mq appears to be linked to 
the central engine. Furthermore, there is a northern molecular gas component related to 
the filamentary dusty spiral seen in the J-H color HST image. This stochastic instability 
might be efficient enough to feed the AGN if it has overcome the likely positive stellar 
gravity torques inside the ILR. 
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